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ABSTRACT 

We study global non-axisymmetric oscillation modes and instabilities in magnetosphere- 
disc systems, as expected in neutron star X-ray binaries and possibly also in accreting 
black hole systems. Our two-dimensional magnetosphere-disc model consists of a Ke- 
plerian disc in contact with an uniformly rotating magnetosphere with low plasma 
density. Two types of global overstable modes exist in such systems, the interface 
modes and the disc inertial-acoustic modes. We examine various physical effects and 
parameters that influence the properties of these oscillation modes, particularly their 
growth rates, including the magnetosphere field configuration, the velocity and density 
contrasts across the magnetosphere-disc interface, the rotation profile (with Newto- 
nian or pseudo General Relativistic potential), the sound speed and magnetic field 
of the disc. The interface modes are driven unstable by Ray leigh- Taylor and Kelvin- 
Helmholtz instabilities, but can be stabilized by the toroidal field (through magnetic 
tension) and disc differential rotation (through finite vorticity). General relativity in- 
creases their growth rates by modifying the disc vorticity outside the magnetosphere 
boundary. The interface modes may also be affected by wave absorption associated 
with corotation resonance in the disc. In the presence of a magnetosphere, the inertial- 
acoustic modes are effectively trapped at the innermost region of the relativistic disc 
just outside the interface. They are driven unstable by wave absorption at the corota- 
tion resonance, but can be stabilized by modest disc magnetic fields. The overstable 
oscillation modes studied in this paper have characteristic properties that make them 
possible candidates for the quasi-periodic oscillations observed in X-ray binaries. 

Key words: accretion, accretion discs - hydrodynamics - MHD - waves - instabilities 
- X-ray: binaries. 



1 INTRODUCTION 

Quasi-periodic oscillations (QPOs) in X-ray fluxes, with frequencies comparable to the dynamical frequencies of neutron stars 
(NSs) and black holes (BHs), have been observed since the 1990s in many X-ray binary systems (e.g., van der Klis 2006), 
largely thanks to NASA's Rossi X-ray Timing Explorer (Swank 1999). They are of great interest because they probe the 
dynamics of the innermost accretion follows around NSs and BHs, and may potentially help constrain the physics of dense 
nuclear matter and strong gravity. 

The phenomenology of kHz QPOs in NS low-mass X-ray binaries is well established (van der Klis 2006). They usually 
occur in pairs, with frequencies between 300 Hz and 1.2 kHz, both varying significantly as a function of the X-ray flux. There is 
some correlation between the upper frequency u u and the lower frequency vi. Spectral analysis indicates most kHz QPOs arise 
from variations in the inner disk boundary layers (Gilfanov et al. 2003). High-frequency QPOs (HFQPOs) (40-450 Hz) in BH 
X-ray binaries have also been intensively studied (Remillard & McClintock 2006; Belloni et al. 2011; Altamirano, Belloni & 
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Figure 1. Schematic description of the cylindrical magnetosphere-disc model considered in this paper. The thick black line on the left 
represents the central compact object. Right next to it is a region of low density plasma threaded by both vertical and toroidal magnetic 
fields. The disc (on the right) consists of fluid of high density (compared to the magnetosphere) and zero or weak magnetic field. These 
two regions are separated by a thin interface/boundary layer. 



Linares et al. 2011; Altamirano & Belloni 2012), although the signals are much weaker and transient. They are only observed 
in the "transitional state" (or "steep power law state") of the X-ray binaries, and have low amplitudes and low coherence. 
Their frequencies do not vary significantly in response to sizeable (factors of 3-4) luminosity changes. Several systems show 
pairs of QPOs with frequency ratios close to 2 : 3 (Abramowicz & Kluzniak 2001), 5 : 3 (Kluzniak & Abramowicz 2002), and 
5 : 2 (Rebusco, Moskalik & Kluzniak 2012). 

Despite much observational progress, the physical origin of HFQPOs remains elusive. Based on the general consensus 
that HFQPOs are associated with the dynamics of innermost accretion flows, a number of ideas/models with various degrees 
of sophistication have been proposed or studied. The general notions of hot blobs exhibiting (test-mass) orbital motion and 
relativistic precession (Stella et al. 1999), nonlinear resonances (Abramowicz & Kluzniak 2001), and in the case of NS systems, 
beats between the orbital frequency and stellar spin (Miller et al. 1998), have been the most popular among the observers, 
although it is not clear how these generic ideas are realized in real fluid dynamical models of accretion flows. Other ideas 
for NS kHz QPOs are discussed in Bachetti et al. (2010). In the case of HFQPOs in BH X-ray binaries, the types of models 
that go beyond the "test-mass motion" idea are those based on global oscillations of accretion flows, including discs, tori 
and boundary layers [see below; see also Sect. 1 of Lai & Tsang (2009) for a critical review]. Systematic studies of accretion 
flow dynamics, combining semi-analytic calculations (for idealized systems to extract physical insights/ingredients) and full 
numerical simulations (for systems with increasing sophistications and realisms) are essential for making further progress in 
this field. 

In this paper we study a class of models involving global oscillations of inner accretion discs and magnetosphere boundary 
layers around NSs or BHs. Our model setup (see Fig. 1) consists of an uniformly rotating magnetosphere with both toroidal and 
poloidal magnetic fields and low plasma density, surrounded by a geometrically thin (Keplerian) disc with high fluid density 
and zero (or low) magnetic field. To make semi-analytical linear perturbation analysis feasible, our model is two-dimensional 
with no vertical dependence (i.e., it is height averaged). We consider both Newtonian discs and general relativistic (GR) discs 
(modelled by the pseudo-Newtonian potential) - we will see that GR can qualitatively affect the linear oscillation modes of the 
system by changing the disc vorticity profile. This kind of magnetosphere-disc systems can form in accreting NS systems when 
the NS magnetic field holds off the accretion disc at a certain distance from the stellar surface (e.g., Ghosh & Lamb 1978). 
They may also be applied to accreting BH systems when magnetic fields advect inwards in the accretion disc and accumulate 
near the inner edge of the disc (e.g. Bisnovatyi-Kogan & Ruzmaikin 1974, 1976; Igumenshchev, Narayan & Abramowicz 2003; 
Rothstein & Lovelace 2008, McKinney, Tchekhovskoy & Blandford 2012). It is usually thought that during the "transitional 
state" (when HFQPOs are observed) of BH X-ray binaries, the accretion flow consists of a thin, thermal disc, truncated by 
a hot, tenuous (and perhaps highly magnetized) inner corona (e.g., Done et al. 2007; Oda et al. 2010). While the true nature 
of this flow is unclear, our simple magnetosphere-disc setup may represent an idealization of such a flow. 

Loosely speaking, two types of global oscillation modes exist in our 2D magnetosphere-disc systems: mertial- acoustic 
modes of the disc and interface modes. Both types of oscillations can be overstable, but driven by different physical mechanisms. 
Our work represents an extension of related previous studies, which have focused on disc oscillations or interface oscillations 
separately. 
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(i) Disc Oscillations: Although relativistic fluid discs can support different types of oscillation modes (Okazaki, Kato & 
Fukue 1987; Nowak & Wagoner 1991; see Wagoner 1999, Kato 2001 and Kato et al. 2008 for reviews), our focus in this paper, 
as well in our previous papers (Lai & Tsang 2009; Tsang & Lai 2008, 2009c; Fu & Lai 2011a), is on the inertial-acoustic modes 
(also called p-modes). These modes have no vertical structure (i.e., their wave functions, such as the pressure perturbation, 
have no node in the vertical direction), and are adequately captured by our 2D model. 

On the other hand, our magnetosphere-disc model does not capture other discoseismic modes, including g-modes (also 
called inertial-gravity modes), c-modes and vertical p-modes (see Wagoner 1999; Kato et al. 2008; Kato 2011a, b). All these 
modes have vertical structures. Non-axisymmetric g-modes and c-modes are either damped due to corotation resonance (Kato 
2003a, Li, Goodman & Narayan 2003; Tsang & Lai 2009a) or their self-trapping zones are easily "destroyed" by disc magnetic 
fields (Fu & Lai 2009). Although axisymmetric g-modes (m = 0) can be resonantly excited by global disc deformations through 
nonlinear effects (Kato 2003b, 2008; Ferreira & Ogilvie 2008), their trapping zones disappear even for a weak (sub-thermal) 
disc magnetic field. By contrast, the basic wave properties (e.g., propagation diagram) of p-modes are not strongly affected 
by disc magnetic fields (Fu & Lai 2009), and it is likely that these p-modes are robust in the presence of disc turbulence (see 
Arras, Blaes & Turner 2006; Reynolds & Miller 2009). 

Our previous papers (Lai & Tsang 2009; Tsang & Lai 2008, 2009c; Fu & Lai 2011a) have shown that non-axisymmetric 
p-modes trapped in the inner-most region of a BH accretion disc are a promising candidate for explaining HFQPOs in BH 
x-ray binaries because their frequencies naturally match the observed values without fine-tuning of the disc properties, and 
because they become overstable due to wave absorption at corotation point when the general relativistic effect on the disc 
rotation profile is taken into account. Although a toroidal disc magnetic field tends to suppress the instability (Fu & Lai 
2011a), a large-scale poloidal field may enhance the instability (Tagger & Pallet 1999; Tagger & Varniere 2006). Most of these 
previous papers adopt a reflective boundary condition at the inner edge of the disc - this is important for the trapping of disc 
p-modes (see Lai & Tsang 2009). How such a reflection can be achieved was not clear. We will show in this paper that the 
disc-magnetosphere boundary serves as a "reflector" for the disc p-modes. 

(ii) Interface Oscillations: Li & Narayan (2004) were the first to consider the interface modes between a magnetosphere 
(with a vertical magnetic field) and an incompressible disc in the cylindrical approximation (i.e., no z-dependence) . They 
showed that the interface modes can be subject to Ray leigh- Taylor instability and/or Kelvin-Helmholtz instability, depending 
on the density contrast and velocity shear across the interface. The mode frequencies roughly scale as mf2i n , where m = 
1, 2, 3 . . . are the azimuthal mode number and Q- m is the angular frequency of the disc flow at the interface, making the 
interface modes a viable candidate for explaining HFQPOs when the interface radius n n is suitably adjusted. Tsang & Lai 
(2009b) generalized the analysis of Li & Narayan (2004) by considering compressible discs. They showed that a relatively 
large disc sound speed is necessary to overcome the stabilizing effect of disc differential rotation and thereby maintain the 
mode growth. Besides linear calculations, numerical simulations of the magnetosphere-disc interface have also been presented 
in a number of papers (see Kulkarni & Romannova 2008; Romanova, Kulkarni & Lovelace 2008 and references therein). 

One ingredient that has been missing from both Li & Narayan (2004) and Tsang & Lai (2009b) is toroidal magnetic 
fields in the magnetosphere, which could be a very important component (Ghosh & Lamb 1978; Ikhsanov & Pustil'nik 
1996). In this paper, we generalize previous calculations of global non-axisymmetric modes confined near the interface of the 
magnetosphere-disc system by taking into account toroidal magnetic fields in the magnetosphere. We aim to examine the 
possibility of large-scale Rayleigh-Taylor/Kelvin-Helmholtz instabilities of the interface in the presence of magnetic fields and 
disc differential rotation. 

Overall, the goal of this paper is to present a complete analysis of large-scale non-axisymmetric modes (including their 
instabilities) in the magnetosphere boundary and in the surrounding disc, and to assess their viabilities for explaining HFQPOs 
in accreting compact binaries. 

We note that in addition to Li & Narayan (2004) and Tsang & Lai (2009b), our study on the interface modes complements 
several other works on the instability of magnetized accretion discs. For example, Lubow & Spruit (1995) and Spruit, Stehle 
& Papaloizou (1995) considered the magnetic interchange instability of a thin rotating disc when a poloidal magnetic field 
provides some radial support in the disc. Lovelace, Turner & Romanova (2009) studied the instability of a magnetopause for 
cases where the shear layer has appreciable radial width and found that the Rossby wave instability may arise in the shear 
layer. Lovelace, Romanova & Newman (2010) considered a setup similar to the present paper, but they included only vertical 
field in the magnetosphere, and focused on small-scale (with radial wavelength <C r) Kelvin-Helmholtz modes in the boundary 
layer. 

This paper is the fourth in our series of papers on the effects of magnetic fields on the global instabilities of various 
astrophysical flows, with the previous three focusing on black hole accretion discs (Fu & Lai 2011a), accretion tori (Fu & Lai 
2011b) and rotating proto-neutron stars (Fu & Lai 2011c), respectively. Our paper is organized as follows. In Section 2, we 
present the setup of our cylindrical magnetosphere-disc model and the dynamical equations for wave modes in the model. 
Section 3 focuses on the interface modes when Newtonian potential is used for disc rotation, while Section 4 deals with both 
the interface modes and disc p-modes when the GR effect (using pseudo-Newtonian potential) is included in the disc rotation. 
We discuss our results and conclude in Section 5. 
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2 EQUILIBRIUM AND PERTURBATION EQUATIONS 

We consider a magnetosphere-disc model similar to the one in Tsang & Lai (2009b). It consists of a magnetosphere region 
where magnetic pressure dominates over gas pressure and a disc region which has high density compared to the magnetosphere 
(see Fig. [TJ. These two regions are separated by an interface (boundary layer). Unlike Tsang & Lai (2009b), who considered 
only vertical magnetic field for the magnetosphere, we take into account both vertical and toroidal fields. In the disc region, 
gas pressure dominates over magnetic pressure. Since any initial poloidal field is likely to generate a dominating toroidal field 
due to the disc differential rotation, for simplicity we take the disc to be threaded by toroidal B field only. We assume that 
flows in both regions are non-self-gravitating, satisfying the usual ideal MHD equations 

^ + V-(pv) = 0, (1) 

^- + {v- V)w = -ivn-V* + -^-(B-V)B, (2) 

at p 4:Tvp 

^- = Vx(vxB), (3) 

where p, P, v are the fluid density, pressure and velocity, $ is the gravitational potential due to the central compact object, 
and 

n.P+f: (4) 

is the total pressure. The magnetic field B also satisfies the equation V ■ B — 0. 

We adopt the cylindrical coordinates (r, <f>, z) which are centered on the central object and have the z-axis in the direction 
perpendicular to the disc plane. The unperturbed background disc has a velocity field v = rf2(r)</>, and the magnetic field 
may consist of both toroidal and vertical components B = B < / > (r)(f> + B z (r)z. The gravitational acceleration in radial direction 
is defined as 

g= ~r (5) 
so that — V3> = —gf and g — rQ^ > 0, where Qk is the angular frequency for a test mass (the Keplerian frequency). Thus 
the radial force balance equation reads 

o2 dn B l t*\ 

pg = pr n---— ( 6 ) 

To investigate the dynamical properties of the flow, we perturb the MHD Eqs. Jf)-© by rewriting any physical variable 
/ as / + Sf with Sf <5cC |/|. Since the unperturbed state is axisymmetric and steady, we consider all perturbation variables 
having the form Sf oc e lm <t>~ lhJt ^ w ith m being the azimuthal mode number and to the wave frequency. Note that the background 
flow and magnetic field have no ^-dependence and we assume that the perturbations also have no z-dependence. The resulting 
linearized perturbation equations contain the variables Sv, Sp, 5P, 5T1 and SB. For mathematical convenience, we define a 
new variable 

5U SP B SB 

Sh= — = — + — . 7 

p p Aitp 

Moreover, using Av — Sv + £ ■ V« = d£/dt = + (v ■ V)£, we find that the Eulerian perturbation Sv is related to the 

Lagrangian displacement vector £ by Sv = — iui^ — rQ. 1 '£ r 4>, where prime denotes radial derivative and 

uj = oj — mQ, (8) 

is the Doppler-shifted wave frequency. In the next two subsections, we will combine the perturbations equations into two 
first-order differential equations (ODEs) in terms of £ r and Sh for the magnetosphere region and the disc region, respectively. 



2.1 The magnetosphere 

In the magnetosphere region (r < r- m ), the flow is assumed to be incompressible and have uniform density (p is constant). For 
this particular magnetized fluid system, the detailed linearized perturbation equations have been given in Fu & Lai (2011b). 
Here we just display the final two ODEs for and Sh: 

^ = Antr + A 12 5h, (9) 
dr 

^ = A 21 i r + A 22 5h, (10) 
dr 

where 



© 2012 RAS, MNRAS 000,ITJfT7l 



Dynamics of the innermost accretion flows around compact objects 5 



An = — 



J A<t> 



At 



A 2 i = . 



2rQ. h 

dr 



.din B. 



dlnr 



- 1 



2 



a 22 



2m wfi + mw 



A<4 



(11) 

(12) 
(13) 

(14) 



and tJA^ = VA<t> /r = B^ / {ry/Anp) is the toroidal Alfven frequency. Note that although we assume low plasma density for the 
magnetosphere, we still require that density is not too low in order to keep Alfven speed far less than the speed of light. 
Otherwise, our original MHD equations break down (e.g., Lovelace, Romanova & Newman 2010). Also note that Eqs. @ and 
(|10f) are the same as Eqs. (119) and (120) (derived for a purely toroidal magnetic field) of §83 in Chandrasekhar (1961). The 
vertical magnetic field B z does not appear in our perturbation equations because we assumed k z = 0, i.e., vertical field lines 
are not perturbed. 



Defining a — Co 



can further combine Eqs. © and (| lOf) into one single equation 



1-rrC 



-ir = 0. 



Two special cases are of interest: 

(i) For systems with B^ = 0, Eq. (|15[) reduces to 



r cj 
or equivalently 



<r = 0, 



r r 
where W = rSv r and 



r d 
mCo dr 



w = o, 



2Q d a 
r dr 



nl/2 



(15) 



(16) 



(17) 



(18) 



is the radial epicyclic frequency. Note that Eq. (|17[1 recovers Eq. (12) in Tsang & Lai (2009b). For either uniform rotation 
profile [fi = const, k = 20., thus (k 2 /2S1)' = 0] or uniform angular momentum profile (Q, oc r -2 , thus k = 0), it has the same 
simple solution 

W oc r m or Sv r oc r m_1 . (19) 
The corresponding solution for £ r is 



t,r oc r 



if Q — const; ^ r oc 



mfi(r) 



if fl oc r 



(20) 



(ii) In the special case of Q = const and B<f, oc r so that ojac/> oc B^jr^fp is constant (note that we have assumed constant 
density in the magnetosphere), [ln(r 3 <7 2 )]' in Eq. (|15[) reduces to [ln(r 3 cD 2 )]'. Therefore magnetic fields in Eq. (|15|) (which 
enters through term a 2 ) completely disappear and the equation is exactly the same as the one for Q. = const and B^ = 0. 
So the solutions for the magnetosphere region in this case are also £ r oc r m ~ , Sh oc r m and the relation between these two 
wavefunctions can be obtained by substituting them back into Eqs. ((9))- (|10[) . Our calculations will focus on this particular 
magnetosphere setup. 



2.2 The disc 



In the disc region (r > n n ), we assume the flow is barotropic so that SP = c? s 5p with c s = \J dP/dp being the sound speed, 
which will be parametrized by a constant c s = c a /(rQ) in our computation. For concreteness, we will also assume a power-law 
disc density profile p oc r~ v with p being a constant. For simplicity, we take the disc toroidal magnetic field to be uniformly 
distributed, i.e., B^ = const. From Eq. ([6]), we obtain the disc rotation profile as 

- -Sftr (2D 



pel 
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In the above expression, we have ignored the effect of magnetic field as the disc field is supposed to be fairly weak in our model 
(see Fig.[T]). The perturbation equations in component form can also be presented in the same form as Eqs. <(9j and (|10[) with 
the detailed expressions of An, A12, A21 and A22 given in Fu & Lai (2011a). For a non- magnetic disc, these equations reduce 
to 



dr 

dSh 
dr 



2mQ dln(rp) 



l ~2 



dr 

2mQ, 



+ 



£r + 



Sh. 



Sh, 



(22) 
(23) 



2.3 The interface 

In the equilibrium state, pressure balance at the interface reads 

Pm + Pbz + Pb*=Pd, (24) 

where P m , Pbz and Pb(p are the magnetosphere gas pressure, magnetic pressure of B z and magnetic pressure of just inside 
the interface, respectively, while Pd is the disc total pressure just outside the interface. This equality imposes an upper limit 
on the strength of the magnetosphere toroidal B field, Pb^/Pd < 1. Defining b — (ujA4>)\r in /^d and p = (pd — p m )/(Pd + Pm), 
where ild is the disc rotation rate at the interface, p m and pd are the fluid densities of the magnetosphere and the disc at the 
interface, respectively, then the above inequality becomes 

P d ~ (1 + p)2c? s 1 °> 

Note that in the magnetosphere, magnetic pressure dominates over gas pressure (Pb z +Pb</> S> Pm) - So Pb<p/Pd — Pb<p/ (Pbct>+Pbz) 
approximately characterizes the relative strength of the toroidal field compared to the vertical field. 

In the perturbed state, we demand that both the radial Lagrangian displacement £ r and the Lagrangian perturbation of 
total pressure All = SH + £ r II' be continuous across the interface. The latter of these requirements yields 

Pm [Sh + £ r r(Q. 2 - Ok - w^)] r . n _ = Pd [Sh - £ r rpc 2 Q. 2 ] . (26) 

In the magnetosphere, we have the analytical solutions for Sh (oc r m ) and £ r (oc r m_1 ) (see Sec. 2.1), which are related via 

Sh = — [u 2 + 2Cjfl + m(2 - rn)uj 2 A4 ] £ r . (27) 

We substitute this relation into the left hand side of Eq. (|26[) and consolidate it with another requirement (continuity of £ r ) 
to obtain the matching condition across the interface 

(w - mf2 m )[w - (m - 2)Q m ] 2 1 + p Sh , [ n 2^pc^ , iNjl2 



+ O z - ! £ _ + 

(1 - p)r in ?r 



+ (m- 1)6 



(28) 



where 5fe and £ r are disc solutions at the interface, Q m and Qd are the rotation rates of the magnetosphere and the disc at 
the interface, respectively. In the case of B<t, — in the magnetosphere, this matching condition reduces to Eq. (22) in Tsang 
& Lai (2009b). Note that Eq. (|28jl has no explicit dependence on the disc magnetic field. This has to do with the fact that 
we chose an uniformly distributed disc B field. 



3 INTERFACE MODES: NEWTONIAN POTENTIAL 

Li & Narayan (2004) and Tsang & Lai (2009b), who considered incompressible and compressible discs respectively, have shown 
that interface modes of magnetosphere-disc boundary can be subject to Ray leigh- Taylor and/or Kevin-Helmholtz instabilities, 
depending on the density contrast and velocity shear across the interface. In order to investigate the effects of magnetosphere 
toroidal magnetic fields on these instabilities, in this section we assume the disc region to be magnetic- field free and employ 
the standard shooting method (Press et al. 1992) to solve Eqs. (122[) - (123[) using the matching condition Eq. ((28} as the inner 
boundary condition at ri n and outgoing wave condition (Tsang & Lai 2009b) at the outer boundary of the disc (r out ). The 
complex mode frequency oj can then be determined as an eigenvalue of the system. Note that in general we also need to 
apply shooting method in the magnetosphere region (see Appendix A) , but for the specific magnetic field profile (and uniform 
rotation) of the magnetosphere we chose in this study, the perturbation equations in the magnetosphere have simple analytic 
solutions (see §2.2) and the effect of magnetosphere on the interface modes is embodied by the interface matching condition 
(i.e, inner boundary condition for disc perturbation equations). 

In our computation, we use p = 1.5, c s — 0.15, and put the outer boundary at r out = 2.85ri n . We vary the strength 
of magnetosphere toroidal magnetic field b (therefore Pb$/Pd\ Eq. [25]) for different sets of (Q m /Qd, pm/ Pd) to see how the 
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Figure 2. The wavefunctions for the unstable interface mode with / pd = 1/99, flm/^d = 1 in a magnetosphere-disc system 

without (left) and with (right) magnetosphere toroidal magnetic field, respectively. The vertical long-dashed lines mark the position of 
the interface which separates magnetosphere region and disc region. The solid and dotted lines represent the real and imaginary parts of 
the wavefunctions, respectively. 



eigenvalue uj will be modified. Throughout this section, we use Newtonian potential $ = —GM/r, so that the Keplerian 
frequency is 57k = (GM/r 3 ) 1 / 2 , and the epicyclic frequency k equals Q. 

3.1 Numerical solutions 

Fig. [2] shows two example wavefunctions for the unstable interface mode. The left panel depicts the case with no toroidal 
magnetic field in the magnetosphere (but keep in mind that there is always finite vertical B field) while in the right panel 
the magnetosphere has a non-zero B^. Since there is no velocity shear across the boundary in both cases (Q. rrl /Q.d = 1), 
the instability is of Rayleigh- Taylor type. We see that the addition of toroidal field has a small effect on the eigenfunction. 
However, as the numbers in the figure indicate, the growth rate of the unstable mode is reduced by more than 50% (from 
0.256f2d to O.llfid) even though the the toroidal magnetic field pressure is only 10% of the disc gas pressure (both evaluated 
at the interface), i.e. even though the toroidal field is much weaker than vertical field in the magnetosphere. 

The suppressing effect of toroidal magnetic field on the instability can be more easily seen in Fig. [3] where we plot the 
real and imaginary parts of the eigen-frequencies for the interface mode as a function of B$. We consider four different sets of 
(Q. m /Q.d, pm/ pd) such that the panels in the same row have the same velocity shear but different density contrast while the 
panels in the same column have the same density contrast but different velocity shear. The modes in the top two panels are 
subject to Rayleigh- Taylor type instability as there are no velocity shear at the interface. The bottom two panels, however, 
have non-zero velocity shear, thus are subject to instabilities of both Rayleigh- Taylor and Kelvin-Helmholtz types. In all 
cases, we observe that the inclusion of toroidal magnetic fields can significantly diminish the growth rate of the unstable 
modes, even completely kill the instability. The critical for absolute shut-down of unstable modes depends on the detailed 
interface parameters (sound speed in the disc, velocity shear, density contrast, etc.). In the case with the most unstable 
modes (bottom-right panel) that we have calculated, the required Pb^/Pd to fully turn off the mode growth is close to one. 
This means that the toroidal magnetic field needs to dominate over the vertical field in the magnetosphere (remember that 
Pb<f>/Pd — Pb<i>/{Pb<p + Pbz))- The mode frequencies (real part of uj), on the other hand, are barely affected by £?0 (at least 
when the growth rate is not very small). By comparing panels in the same column, we see that increasing velocity shear at 
the interface leads to larger growth rates, which is not surprising as more velocity shear means more free energy that the 
system can tap on to drive the instability. Comparing panels in the same row reveals an interesting features. We found that 
larger density contrasts leads to smaller growth rates, which apparently contradicts the standard Rayleigh- Taylor instability 
result. In particular, the m = 1 mode is totally stable when pm/pd = 1/99 while becomes unstable when p m /pd = 1/9 and 
B^ does not affect the m = 1 mode at all. Comparing different lines in each individual panel shows that modes with larger m 
generally have larger growth rates, but they are also more easily stabilized by the toroidal field. 
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Figure 3. Real and imaginary parts of the wave frequencies (in units of fid, the disc rotation rate at the interface) for unstable interface 
modes as a function of magnetosphere toroidal magnetic field strength. The top and bottom panels are for the cases without and with 
velocity shear at the interface, respectively, while the left and right panels depict models with the magnetosphere to disc density ratio 
being 1/99 and 1/9, respectively. Different line types are associated with different values of m, with the solid lines for m = 1, dotted 
lines for m = 2, short-dashed lines for m = 3, long-dashed lines for m = 4, and dot-dashed lines for m = 5. 



3.2 Discussion of results 

The zero effect of on the m = 1 mode (whether it is stable or unstable) can be easily understood by looking at the 
matching condition Eq. (|28|l . This equation (i.e., inner boundary condition for disc equations) is the only place that affects 
the eigenvalue problem, and the magnetic effect comes in via the last term on the right hand side of the equation [(m—l)b 2 Q^]. 
When m — 1, this term vanishes so that the equation is the same as in the case with zero B$. 

To explain the other features observed in Fig. [3] we carry out a local analysis of the effects of magnetic fields on Rayleigh- 
Taylor/Kelvin-Helmholtz instability in a plane-parallel flow (see Appendix B). From the last term in the final expression of 
wave frequency (Eq. |Blj ). we see that magnetic tension (k ■ B) provides a suppressing force against Rayleigh-Taylor/Kelvin- 
Helmholtz instability of a two-layered fluid system. This is because extra work needs to be done in order to increase the 
boundary layer deformation (i.e., growing perturbation). Our cylindrical magnetosphere-disc model resembles the simple 
parallel two-layered flow in that our wavenumber in azimuthal direction m/r acts like k in the x direction and our B^ also 
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lies along the direction of the wave vector. Thus, the same stabilizing mechanism also applies in our system, i.e., the tension 
force of the toroidal magnetic field suppresses the magnetosphere-disc interface instability. 

Qualitatively, the growth rate of interface modes in our system is determined by combined effects of four factors: density 
contrast (pm/pd), velocity shear across the interface (Q m /Qa), degree of differential rotation in the disc, and toroidal field of 
the magnetosphere. The first two tend to enhance instability, while the latter two tend to suppress instability. In different 
limiting cases, approximate analytic expressions for the mode growth rate can be derived: Li & Narayan (2004) considered 
a case where the disc is incompressible with constant density and came up with their Eq. (45), while Tsang & Lai (2009b) 
studied a compressible disc whose density is much larger than the magnetosphere density, and summarized the results in their 
Eqs. (26) and (28). Although for the generic conditions studied in this paper, an analytical expression cannot be rigorously 
derived, we can write the mode growth rate schematically as follows: 



IT + W KH - W vort - ^6 ( 29 ) 

with 

Wrt - 2(1 + n)mVt 2 eSA - 2(1 - n)mftl Stm _, (30) 
cjkh — (1 - p 2 )m 2 (Q d - Q, m ) 2 + 2(1 - p 2 )m(Q, d - Q m ){(d - Cm), (31) 

^vort - (Cd - Cm) 2 + 2/i(Cd - Cm) + V (Cd + Cm) 2 , (32) 

m 2 B 2 

W? * ^r\ Pi lp m ) * (1 - ^™ 2 -V2, (33) 

where rQ 2 H = d$/dr — rQ. 2 is the effective gravitational acceleration, C = /t 2 /(2fi) is the fluid vorticity and as before the 
subscripts d and m denote disc and magnetosphere, respectively. Several features can be noted in Eq. Q29p : 

(i) The first two terms under the square root are the Rayleigh- Taylor and Kelvin-Helmholtz types of destabilizing 
factors while the last two terms characterize the stabilizing effects due to flow vorticity and toroidal magnetic field in the 
magnetosphere, respectively. The suppressing effect of finite fluid vorticity has been discussed in both Li & Narayan (2004) 
and Tsang & Lai (2009b). We see that the various terms have different dependences on the density contrast p. As p increases 
(fj, = when p m /pd = 1, and p = 1 when Prn /p d = 0), increases, making the system more Rayleigh- Taylor unstable; 
at the same time, w^h decreases, making the system less Kelvin-Helmholtz unstable. In addition, cavort a l so becomes larger, 
leading to stronger suppressing effect. This explains why in Fig. [3] the mode with p m / Pd = 1/99 is less unstable than the one 
with p m /pd = 1/9. 

(ii) Since Q. 2 S d = (l/r)(d<E>/dr) — Q. 2 = — (c 2 /p)(dp/dr), when c s is too small (i.e., small effective gravity in the disc), the 
destabilizing terms would not be able to compete with the stabilizing terms, resulting in a stable system. Thus, a sufficiently 
high disc sound speed is needed to attain the interface instability (see Tsang & Lai 2009b). 

(iii) When there is no B^, we see that uj 2 ^ an d w|h both depend on m while w 2 ort does not. Thus, modes with higher 
m tend to be more unstable, as shown in Fig. [3] (the Pb<t>/P d = case; see also Li & Narayan 2004 and Tsang & Lai 2009b). 
The m = 1 mode would be even less unstable if the density contrast is too big [see (i) above]. This explains why in Fig. [3] the 
m = 1 mode with p m /Pd = 1/99 is stable while the one with p m / P d = 1/9 is unstable. 

(iv) When 7^ 0, besides cj 2 , the magnetic field strength also appears in where Q^a m contains a term that 
is proportional to ~ B 2 , / Air p m r 2 (note the minus sign). This shows that the toroidal field in the magnetosphere plays two 
different roles in determining the stability of system. On the one hand, the magnetic tension resists perturbation growth, 
thus suppressing any instability; on the other hand, the magnetic force increases the effective gravity (pointing towards the 
center) of the background flow, thus promoting Rayleigh- Taylor type instability. The former effect is proportional to m while 
the latter to m 2 . Hence, in general the suppressing effect is more important. This is consistent with the fact that in Fig. [3] 
the growth rates of modes with higher m decrease more rapidly as B^ increases. As noted above, B<f, does not affect the 
m = 1 mode because the last term in Eq. (|28[) vanishes for m = 1. Now we have a better understanding of what this means 
physically: For the m = 1 mode, the aforementioned two opposing effects associated with B^, happen to cancel each other. 
This exact cancellation, however, cannot be captured by the approximate expression Eq. (|29p . 

Overall, we see that Eq. Q29p . though schematic, is quite useful in explaining most of the numerical results of this section. 
Note that various flow parameters, such as the density contrast p, azimuthal mode number m, toroidal magnetic field B^ 
and disc vorticity Cd, appear in more than one of the four terms in this equation. Thus, they are associated with both the 
stabilizing and destabilizing effects. The only exception is perhaps the velocity shear at the interface, which always facilitates 
the interface instability. 

Finally, we note that the interface instability associated with magnetosphere-disc boundary studied in this paper is 
qualitatively different from those of Lubow & Spruit (1995) and Spruit, Stehle & Papaloizou (1995), who carried out local 
analysis of a thin rotating disc threaded by a large-scale poloidal field. Lovelace, Romanova & Newman (2010) considered 
similar interface instability as in our paper, but they focused on small-scale modes and included only vertical field in the 
magnetosphere. 
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4 INTERFACE MODES AND DISCO-SEISMIC MODES: PSEUDO-NEWTONIAN POTENTIAL 

In this section, we consider how the global oscillation modes in our disc-magnetosphere system (Fig. 1) are modified by general 
relativistic (GR) effects. In particular, GR changes the disc rotation profile Q.k(j) and makes the radial epicyclic frequency 
k(t) a non-monotonic function of r: As r decreases, k first increases, attains a maximum value and then falls to zero at the 
Innermost Stable Circular Orbit (ISCO), nsco = 6GM/c 2 (for a non-spinning compact object). If the inner disc boundary is 
close to nsco, then this non-monotonic k profile can have two consequences: (1) It significantly reduces disc vorticity (k 2 /2S1) 
near the inner disc boundary, therefore helps the unstable interface modes grow as there are less disc vorticity suppressing 
effects (see the discussion in Section 3.2; see also Tsang & Lai 2009b); (2) It could lead to a vortensity (k 2 /2Q,p a.k.a potential 
vorticity) profile that has positive gradient in the inner disc region. As found in our previous studies (Tsang & Lai 2008; Lai 
& Tsang 2009; Fu & Lai 2011a), positive vortensity gradient renders the disc inertial-acoustic modes (p- modes) unstable. In 
Section 3, we chose a density profile p oc r~ 3 ^ 2 so that the vortensity profile (with the Newtonian potential) is completely flat, 
and thus we only found unstable (overstable) interface modes. In this section, with the GR effect included, both the interface 
modes and disc p-modes can be unstable and they co-exist in our disc-magnetosphere system. The main goal of this section 
is then to study these two modes and the effects of the magnetosphere toroidal B field, disc toroidal B field and inner disc 
boundary location. 

We employ the pseudo-Newtonian potential (Paczynski & Wiita 1980) to mimic the GR effect: 

* = --^-, (34) 

r - r s 

where rs = 2GM/c 2 is the Schwarzschild radius. The corresponding Keplerian rotation frequency and epicyclic frequency are 



n K = (35) 



K = nJ- — — . (36) 

V r — rs 

Except for Section 4.3, we will take the magnetosphere-disc interface to be located at the ISCO, n n = nsco = 3rs. 



4.1 Non-magnetic discs 

We first consider the case where the disc outside the magnetosphere is non-magnetic. We solve the same equations and use the 
same interface boundary condition as in Section 3 except that we replace the Newtonian potential with Pseudo-Newtonian 
potential. 

Figure [4] shows the complex eigenfrequencies of disc inertial-acoustic modes and interface modes as a function of the 
dimensionless magnetosphere toroidal field strength Pb^/Pd- The other parameters are indicated in the figure. As noted 
above, because of the GR effect, the disc vortensity profile has a positive slope near the ISCO (see Fig. [5}. The disc inertial- 
acoustic mode (p-mode) of lowest radial order (for a given m) has Re(o;) < mQ.d, and the corotation resonance lies in the 
disc (see the right panels of Fig. [6] for an example). Corotational wave absorption then makes these modes unstable. We see 
from Fig. [4] that as long as there is a magnetosphere to serve as an inner boundary condition for the disc mode, the mode 
frequency and growth rate depend rather weakly on Pb^/Pd- 

The interface mode is also strongly influenced by the the GR effect. Comparing the right panels of Fig. [4] with the 
bottom-left panel of Fig. [3] we see that for Pbcf, = 0, the mode growth rate (for a given m) is larger when the GR effect is 
included. Again, this is because in GR, the disc vorticity is smaller near the ISCO, leading to less rotational suppression of the 
RT instability (see Section 3.2). As Pb<f>/Pd increases, the interface mode growth rate first decreases (just as in the Newtonian 
case) due to magnetic tension, then starts to increase beyond certain critical value of Pb<f,/Pd (for a given m). This behaviour 
arises from the effect of corotation resonance (see Fig. [5}: For small Pb<f>/Pd, the interface mode has Re(aj) > mfld, thus 
no corotation resonance exists in the flow; as Pb^/Pd becomes larger, the mode frequency drops below mfld, and corotation 
resonance comes into play, which overwhelms the suppression effect of the magnetic tension (see the left and middle panels 
of Fig. [6] for two examples of mode wavefunctions) . The growth rate of such interface mode (labelled as Type II in Fig. [5)l 
is much larger than the corresponding disc p-mode because the corotation resonance lies much closer to r- ln for the interface 
mode than for the p-mode - this gives rise to much stronger wave absorption at the corotation resonance. 



4.2 Magnetic discs 

Now we consider the case in which the disc outside the magnetosphere has a finite toroidal magnetic field. For simplicity, we 
take the limit that the density in the magnetosphere p m goes to zero. In this case, the interface boundary condition Eq. (|28|l 
reduces to ALT = 0, i.e., the Lagrangian perturbation of total disc pressure equals zero at Ti n . With this boundary condition, 
the setup is the same as in Fu & Lai (2011a) and we only need to consider the dynamical equations for the disc. 
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Figure 4. Real and imaginary parts of the wave frequencies (in units of Cl d , the disc rotation rate at rj n ) for unstable inertial-acoustic 
disc modes (left panels) and unstable interface modes (right panels) as a function of the magnetosphere toroidal field strength. Different 
line types represent different azimuthal mode number m. The sound speed profile, density contrast and velocity shear across the interface 
are the same as in the bottom- left panel of Fig. [3] The only difference from Fig. [3] is that GR effect (using pseudo-Newtonian potential) 
is included in the disc rotation profile, which renders both types of modes unstable. On the right panels, the vertical lines (for m = 2 
and 3) mark the point when corotation resonance moves into the flow (see Fig. 0. 



Again, two types of oscillation modes exist in our system (see Fig. [Jj: the disc inertial-acoustic modes and interface 
modes. The effect of the disc magnetic field on the inertial-acoustic modes has already been studied in detail by Fu & Lai 
(2011a). We see from the left panels of Fig. [Jjthat the mode frequency is only slightly modified by the disc B^, but the growth 
rate is reduced so that the mode become stable even for modest (sub-thermal) disc toroidal fields. 

The interface mode (the right panels of Fig. [7]) was not considered by Fu & Lai (2011a). We see that as the disc magnetic 
field increases, the mode frequency varies modestly (about 10%) for a wide range of B^, but the growth rate changes more 
significantly. The dependence of the growth rare as a function of disc B^ can be understood as follows: (i) For B^ — 0, 
RT instability drives the mode growth; (ii) As B$ increases, the magnetic tension tends to suppress the growth; (iii) As B^, 
increases, the real mode frequency decreases. When B$ exceeds some critical value, corotation resonance appears in the flow 
(disc), and wave absorption at corotation then enhances the growth rate. But as B^ increases further, the corotation radius 
lies at a larger distance from rj n , thus the corotational effect becomes less important (because of the large evanescent zone 
separating r- m and the corotation radius) and the mode growth rate decreases again also because of magnetic tension. 

Concerning (iii) above, as noted in Fu & Lai (2011a), in the presence of the disc toroidal magnetic field, the corotation 
resonance (where Cj = 0) is split into the inner/outer magnetic resonances, where 

uj = ±mtOA^, (37) 

where loa^, = VA<p/r = B^/(r^/47rp) is the toroidal Alfven frequency of the disc. When these magnetic resonances exist in the 
flow (disc) , wave absorption comes into play in the mode growth. Note that the signs of wave absorption at the two magnetic 
resonances are different. As Re(aj) decreases (with increasing -B^), the outer magnetic resonance (where u = mu)A<i>) first 
enters the flow, causing the mode growth rate to increase. When the inner magnetic resonance (where to — —mujA<t>) enters 
the flow, the wave absorptions at the two magnetic resonances partially cancel, and the mode growth rate starts to decrease 
again (see the right panels of Fig.[7|). In this case, the term "corotation resonance" simply refers to the combined effect of two 
magnetic resonances. 



4.3 Effects of different inner disc radii 

Here we consider the same setup as in Sections 4.1-4.2 except that the inner disc boundary lies outside the ISCO. The 
motivation for considering n n > nsco is that in real accreting NS systems, the magnetosphere radius may well be outside 
the ISCO (e.g., in accreting millisecond X-ray pulsars with surface magnetic field of 10 8-9 G, the Alfven radius is about 1.5-2 
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Figure 5. Wave propagation diagram (upper panel) for both interface modes and disc inertial-acoustic modes accompanied by the disc 
vortensity profile (lower panel). In the upper panel, the three thick solid curves depict the disc rotation profile f2 and f2 ± n/m, where m 
is the azimuthal mode number and ft is the radial epicyclic frequency. Note that since ft(nsco) = ® these three curves join each other at 
the Innermost Stable Circular Orbit (which is also the disc inner boundary r; n in our setup of Sections 4.1 & 4.2). The dashed horizontal 
line marks the point when the corotation radius [where Re(u)/m = Q] happens to be exactly at the inner disc boundary. For modes with 
Re(u)/m above this line, there is no corotation in the flow; for any modes with Re(cj)/m below this line, the corotation exists in the flow 
(fc > r- ln ). Three types of wave modes (labeled as I, II and III) are shown in this diagram where the horizontal solid lines indicate the 
evanescent zones and the wavy curves denote the wave propagation zones. Note that for both Type I and II waves, there is also a tiny 
wave zone (hardly visible in the figure) in the disc region just outside the ISCO. Type I and II waves are both interface modes while 
type III waves are disc inertial-acoustic modes. Note that the magnctosphcrc region (r < rjgco) an d the- region between the inner and 
outer Lindblad resonances [where Re(u;)/m = Q ± K/m] are wave evanescent zones. 




Figure 6. Wave functions for the three types of wave modes depicted in Fig. [5] the interface mode with no corotation resonance in the 
flow (Type I, left panels), interface mode with corotation resonance in the flow (Type II, middle panels) and disc inertial-acoustic mode 
(Type III, right panels). The azimuthal mode number is m = 3 and all other parameters are the same in Fig. [4] The solid and dotted 
lines show the real and imaginary parts, respectively. In the middle and right panels, the vertical dot-dashed lines represent the location 
of corotation resonance. 
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Figure 7. Real and imaginary parts of the wave frequencies (in units of Cl d , the disc rotation rate at rj n ) for unstable inertial-acoustic 
disc modes (left panels) and unstable interface modes (right panels) as a function of the disc toroidal field strength. The x-axis specifies 
the ratio of the Alfven velocity to disc rotation velocity at the inner boundary r; n . The magnetosphere density is set to zero. The disc 
parameters are r; n = nscoi m = 2, c s = O.lrfi and p oc r~ . For the interface mode (right panels), the two dotted vertical lines from 
left to right mark the entering of the outer magnetic resonance and inner magnetic resonance into the flow, respectively. 



stellar radii). Also, in the "transitional state" (when HFQPOs are observed) of BH X-ray binaries, the thin thermal disc may 
be truncated at a radius slightly larger than the ISCO (e.g., Done et al. 2007; Oda et al. 2010). 

As in Section 4.2, we assume that the magnetosphere has a negligible density. Figure [8] shows an example (for a non- 
magnetic disc) of how the (complex) mode frequencies depend on r ln . When expressed in units of fid, the real frequencies 
of both interface and inertial-acoustic modes are only modestly affected by n n /risco, but the growth rates decrease rapidly 
with increasing ri n /risco- For the inertial-acoustic mode, this arises from the reduced vortensity slope at distances further 
away from the ISCO, which results in smaller wave absorption at corotation. For the interface mode, the larger vorticity just 
beyond r ln leads to a stronger rotational suppression of the RT instability, thus a smaller mode growth rate. 

Figure [9] shows how the complex frequency of the interface mode depends on the disc magnetic field when the inner 
disc radius is set to 1.5risco- The behaviour of the mode growth rate as a function of (vA<f>/rQ)in is similar to the right panels 
of Fig. [7] Here, the mode growth rate is small when B$ — 0. So when the outer magnetic resonance enters the flow, wave 
absorption dramatically increases the mode growth rate. 



5 DISCUSSION AND CONCLUSION 

In this paper we have studied the non-axisymmetric MHD modes and instabilities in a 2D model of magnetosphere-disc 
systems (see Fig. 1), as may be realized in accreting neutron star or black- hole X-ray binaries (see Section 1). We have 
examined various physical effects and parameters that can influence the global modes in the system, including the density 
and magnetic field of the magnetosphere, the velocity contrast across the magnetosphere-disc interface, the rotation profile 
(Newtonian vs GR), the temperature and magnetic field of the disc. We restrict to modes that do not have vertical structure, 
but otherwise our calculations include all possible instabilities and global oscillations associated with the interface and the 
disc. We highlight several key findings and implications of this paper below. 



5.1 Interface instabilities in a rotating, magnetized system 

This paper includes a comprehensive study of the large-scale Ray leigh- Taylor (RT) and Kelvin-Helmholtz (KH) instabilities 
associated with the interface of a rotating, magnetized system. RT and KH instabilities have been studied intensively in 
plane-parallel flows through both theoretical analysis and laboratory experiments (e.g. Chandrasekhar 1961; Drazin & Reid 
1981), and have found applications in various astrophysical and space environments. But few papers have focused on rotating 
systems (e.g., Spruit et al. 1995; Lovelace et al. 2009, 2010). Our study generalizes previous works by Li & Narayan (2004) 
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Figure 8. The real and imaginary parts of the wave frequencies (in units of Q d , the disc rotation rate at r; n ) for unstable inertial- 
acoustic modes (left panels) and unstable interface modes (right panels) as a function of the inner disc radius rj n (in units of rjgco)- 
The magnetosphere inside rj n has zero density and the disc toroidal magnetic field is set to zero. The other disc parameters are m = 2, 
c s = O.lrQ and p oc r . 
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Figure 9. Real and imaginary parts of the wave frequencies (in units of Q^, the disc rotation rate at ri n ) for unstable interface modes 
as a function of the disc toroidal field strength, with the inner disc boundary located at r; n = 1.5rigco- The rr-axis specifies the ratio 
of the Alfven velocity to disc rotation velocity at r- ln . The magnetosphere density is set to zero. The disc parameters are rj n = nsco> 
m = 2, c a = 0.1rf2 and p oc r~ 1 . The two dotted vertical lines from left to right mark the entering of outer magnetic resonance and inner 
magnetic resonance into the flow, respectively. 
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and Tsang & Lai (2009b) by considering compressible fluid, magnetic field and rotation profile of the disc, as well as generic 
field (poloidal and toroidal) configuration of the magnetosphere. 

As in plane-parallel flows, the interface modes are mainly driven unstable by the RT and KH instabilities. Toroidal 
magnetic field tends to suppress the instabilities through magnetic tension. The magnetic field also indirectly affects the RT 
instability by modifying the effective gravity. Except for the m = 1 mode, for which the two opposite effects of cancel (as 
conjectured by Li & Narayan 2004), we find that increasing B^, generally tends to reduce the growth rates of the interface 
modes. Differential rotation (with finite vorticity) also tends to suppress the interface instability. To overcome this suppression 
effect, the disc must have sufficiently large temperature (sound speed) (see Tsang & Lai 2009b). General relativity (GR) can 
significantly affect the growth rate of interface modes because the disc rotation near the ISCO has smaller vorticity in GR 
than in Newtonian theory. 

Another qualitatively new finding of this paper is that corotation resonance can significantly influence the interface 
instabilities. As the toroidal field strength in the magnetosphere or in the disc increases, the real frequency of the interface 
mode falls below mQd (where fid is the disc rotation rate at the interface), and corotation resonance (or its generalization 
to magnetic resonances) appears in the disc. Wave absorption at corotation can then significantly change the interface mode 
growth rate (see Figs. 2] [JJ . 

5.2 Inertial- Acoustic Modes of relativistic Disc 

Our model system (Fig. 1) also accommodates inertial-acoustic modes (p-modes) of relativistic discs. These modes are driven 
unstable primarily by wave absorption at corotation resonance. The magnetosphere-disc interface naturally serves as the inner 
boundary for the disc. Our study in this paper complements our previous works (Lai & Tsang 2009; Tsang & Lai 2009c; Fu & 
Lai 20f fa) by properly treating the inner boundary condition for disc oscillations. Our result shows that the magnetosphere 
behaves as a robust "reflector" for spiral waves in the disc: The p-mode frequency and growth rate do not depend sensitively 
on the property (density, magnetic field) of the magnetosphere. In agreement with Fu & Lai (20f fa), we find that a modest 
disc toroidal field tends to reduce the growth rate of disc p-modes. 

5.3 Implications for High-Frequency QPOs in X-ray Binaries 

The large-scale (m = 1, 2, 3, . . .), overstable oscillation modes studied in this paper may provide an explanation for the high- 
frequency QPOs observed in NS and BH X-ray binaries (see Section f). Obviously, the simplicity of our model setup precludes 
detailed comparison with the phenomenology of QPOs. But we note the following relevant features of disc inertial-acoustic 
modes and interface modes. 

(1) The disc inertial-acoustic modes have frequencies uj — PamQd (with fid the disc rotation rate at n n ), with /3d < f 
(typically ~ 0.5) depending on model parameters (such as disc sound speed) and boundary conditions. If the inner disk 
is located at the ISCO, the mode frequencies can be computed ab initio, and the results are generally consistent with the 
observations of the HFQPOs in black-hole X-ray binaries (see Lai & Tsang 2009; Tsang & Lai 2009c). The problem with these 
modes is that even a weak (sub-thermal) disc toroidal magnetic field can suppress their instabilities (see Section 4 and Fu 
& Lai 201 fa). Although large-scale poloidal fields can enhance the instability under certain conditions (see Tagger & Pallet 
1999; Tagger & Varniere 2006), it is not clear at this point which effects (disc toroidal field vs large-scale poloidal field) will 
dominate. 

(2) The interface modes have frequencies uj — /3imfid, with /3, ~ f . If r- m ~ nsco, the implied QPO frequencies would 
be too high compared to observations. Of course, as discussed before (see Sections I and 4.3), it is certainly possible that 
rin is somewhat larger than nsco in real systems, but then we would lose the predictive power of our calculations (since in 
our model rj n is a free parameter). On the other hand, the interface modes are robustly unstable, driven by the RT and KH 
instabilities, and by the corotation effect, especially when the GR effect is included. In most cases, the growth rates of the 
interface modes are much larger than the disc modes. 

Overall, the results of this paper suggest that if the real accreting NS or BH systems can be approximated by our 
magnetosphere-disc model, the interface oscillations are more likely than inertial-acoustic oscillations to provide an explanation 
for the observed QPOs. This is reasonable for NS systems, and is consistent with spectral analysis of kHz QPOs in NS X-ray 
binaries (e.g., Gilfanov et al. 2003). For BH systems, this would require that the inner disc radius in the "transitional state" 
to be slightly larger than nsco (e.g., Done et al. 2007; Oda et al. 2010). 

Nevertheless, it should be noted that our calculations of global disc modes are based on a fairly primitive model. Several 
potentially important physical ingredients have been neglected in our setup, such as finite vertical wavenumber, disc turbulence 
due to MRI, radial mass inflow across ISCO, etc.. The effects of these factors on our results are certainly worth future 
investigations. Plus, the disc modes in our study are all linear modes which in principle become nonlinear eventually. Whether 
our results here are still valid during the nonlinear stage or not needs to be addressed by doing numerical simulations (Fu & 
Lai 2012, in prep.). Moreover, we have not done any analysis regarding to how these global modes would manifest themselves 
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observationally. In this sense, our study should better be treated as a demonstration of interesting physical mechanisms, rather 
than a mature theoretical tool for explaining or predicting HFQPO observations. 
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APPENDIX A: GENERAL MAGNETIC FIELD PROFILES IN THE MAGNETOSPHERE 

In the main text of this paper we adopt a special toroidal magnetic field profile B^ oc r in the magnetosphere so that the 
perturbation equations have simple analytical solutions. For a more general magnetic field profile B^, oc r q , with q > 1, 
analytical solutions are no longer attainable, and we must solve numerically the perturbation equations for both the disc 
(Eqs. [22]-[23]) and the magnetosphere (Eqs. [9]-[l0]). To this end, we first need to derive the boundary condition at the 
the inner boundary of the magnetosphere, which is close to the center of the system (r = 0). This is done by requiring the 
solutions to be regular at small radius. As r — > 0, we observe that ou</> oc B^/r — > 0. In this limit, Eqs. ([§))- ()10p then reduce to 

d£ r 1 OJ — 2m£l , m 2 . , „ . 

-J- = = tr + — 5h, Al 

dr r to r 1 

— — = (cj - rQ, 2 Y r H Sh, (A2) 

dr ruj 

i.e, recovering the hydrodynamic equation^]. These simplified equations can be readily solved, leading to a relation between 
£ r and Sh: 

& = rw(£D + 2fi) " (A3) 

With this regularity condition implemented at some finite yet small inner boundary r c , we carry out the integration in the 
magnetosphere towards the interface, while at the same time we also integrate the disc equations towards the same interface 
and employ shooting procedure to solve for the eigenvalues by requiring solutions from both regions satisfy the matching 
condition at the interface (i.e., continuity of £ r and All). Note that the eigenvalues of the system now have two components: 
the wave frequency w, and the relation between the solutions in two regions, e.g., Sh(r c )/5h(r out ) or £ r (r c ) /fr(V out)- 



APPENDIX B: PLANE-PARALLEL FLOW WITH A COMPRESSIBLE UPPER LAYER AND A 
MAGNETIZED LOWER LAYER 

Consider a system with two separate fluid layers in a constant gravitational field g = — gz. The upper layer (z > 0) is a 
compressible fluid of density p = p+e~ z ^ Hz and constant horizontal velocity u+, and the lower layer (z < 0) is an incompressible 
fluid of constant density p_ and horizontal velocity u_. Here, H z — (? s /g is the density scale height with c s being the sound 
speed, and both velocities are along the z-axis. This system was studied by Tsang & Lai (2009b), but here we we add a 
uniform horizontal magnetic field B x in the lower fluid. As in Tsang & Lai (2009b) (correcting typos in their Appendix), 
we apply perturbations of the form e lkx ~ 1 " t to both layers and solve for to by demanding the Lagrangian displacement and 
Lagrangian total pressure perturbation be continuous across the interface between upper and lower fluids. Denoting 

— to — ku+ , 

kl = (k 2 -Hj 2 + /c s ), 



See Fu & Lai (2011c) for the derivation of the regularity condition for a even more general system (differentially rotating, magnetized, 
self-gravitating, etc.). 
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and 



P+ = P+k/k, 



the final solution for uj can be written as 




+ 



Mp+ + p-)' 



k 2 Bl 



(Bl) 



The three terms under the square root clearly have different physics meanings. The first term signifies the Kelvin-Helmholtz 
instability which disappears when velocity shear across the interface vanishes; the second term depicts the Rayleigh- Taylor 
instability which only exists for non-zero density contrast; the last term, the only difference between Eq. (IB1|) and Eq. (A14) 
in Tsang & Lai (2009b), describes the stabilizing effect of the magnetic field that lies along the direction of the perturbation 
wave vector. If we take the incompressible and hydrodynamic limit (i.e., c s — > oo, H z — > oo so that k — > k, p~+ — > p+, and 
B x — > 0), then the above equation reduces to Eq. (B12) in Li & Narayan (2004). 
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